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The low energy neutrino factory has been proposed as a very sensitive setup for future searches 
for CP violation and matter effects. Here we study how its performance is affected when the 
experimental specifications of the setup are varied. Most notably, we have considered the addition 
of the 'platinum' — > channel. We find that, whilst theoretically the extra channel provides 
very useful complementary information and helps to lift degeneracies, its practical usefulness is lost 
when considering realistic background levels. Conversely, an increase in statistics in the 'golden' 
v e — s> Vy, channel and, to some extent, an improvement in the energy resolution, lead to an important 
increase in the performance of the facility, given the rich energy dependence of the 'golden' channel 
at these energies. We show that a low energy neutrino factory with a baseline of 1300 km, muon 
energy of 4.5 GeV, and either a 20 kton totally active scintillating detector or 100 kton liquid argon 
detector, can have outstanding sensitivity to the neutrino oscillation parameters #13, 5 and the mass 
hierarchy. For our estimated exposure of 2.8 x 10 23 kton x decays per muon polarity, the low energy 
neutrino factory has sensitivity to #13 and S for sin 2 (2613) > 10~ 4 and to the mass hierarchy for 
sin 2 (26»i 3 ) > 10- 3 . 

PACS numbers: 14.60.Pq 



I. INTRODUCTION 



Neutrino oscillations have been robustly established. The present data require two large (612 and #23) angles and one 
small (#13) angle in the neutrino mixing matrix, and at least two mass squared differences, Am^ = mf — to 2 (where 
mi's are the neutrino masses), one driving the atmospheric (Amy and the other one the solar (Am^) neutrino 
oscillations. The mixing angles 612 and O23 control the solar and the atmospheric neutrino oscillations, while #13 is 
the angle which connects the atmospheric and solar neutrino regimes. 

A global fit performed at the end of 2008 [l[ (see also @) provides the following 3cr allowed ranges for the at- 
mospheric mixing parameters: jAm^l = (2.07 — 2.75) x 10~ 3 eV 2 and 0.36 < sin 2 623 < 0.67. The sign of Am^ 
(sign(Am| 1 )) cannot be determined from the existing data. The two possibilities, Am|j > or Am|j < 0, correspond 
to two different types of neutrino mass ordering: normal hierarchy and inverted hierarchy respectively. In addition, 
information on the octant of 623, if sin 2 2823 ^ 1, is beyond the reach of present experiments. The best fit values for 
the solar neutrino oscillation parameters are Amfx = 7.65 x 10~ 5 eV 2 and sin 2 6*12 = 0.30 [1]. A non-zero value of #13 
is crucial to enable a measurement of the CP violating phase 8 and the mass hierarchy. A combined three-neutrino 
oscillation analysis of the solar, atmospheric, reactor and long-baseline neutrino data 1] constrains the third mixing 
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angle to be sin 2 #13 < 0.056 at the 3<r confidence level, with a best fit value of 0.01. Different analyses undertaken in 
2008 using the available solar data from the third phase of the Sudbury Neutrino Observatory (SNO-III) and recent 
data from KamLAND favour a non-zero value of sin 2 6*13 at ~ la. Similarly, Super-Kamiokande data on atmospheric 
neutrinos [H leads to a ~ 2a preference for sin 2 6*13 > 0. This second claim is, however, controversial Q. 2009 
data from the MINOS experiment, studying the appearance channel — > u e , also favours a non-zero value of #13 
but with an even larger best fit, even more in conflict with the stringent upper bound that comes mainly from the 
CHOOZ reactor experiment of sin 2 #13 < 0.056. A preliminary combination of all the data provides a la range of 
sin 2 6*i3 = 0.02 ± 0.01 Q. This hint for non-zero #13 and the resulting tension among the different datasets will be 
probed by the forthcoming generation of accelerator [f| @] and reactor [8l-flr3| experiments. If the hint for large #13 
is confirmed, the exciting possibility to search for leptonic CP violation, encoded in the phase 8, and the ordering of 
neutrino masses [Tl| . will be open. However, these experiments lack the required sensitivity and a new generation 
of neutrino oscillation experiments is therefore needed for this task, and to explore even smaller values of #13 if the 
present hint is not confirmed. 

Future long-baseline experiments will require powerful machines and extremely pure neutrino beams. Among these, 
neutrino factories (l2^ . in which a neutrino beam is generated from muons decaying within the straight sections of 
a storage ring, have been shown to be sensitive tools for studying neutrino oscillation physics [l2| - [27| . The neutrino 
factory exploits the golden signature of the wrong-sign muon (l2l [l3T | events, i.e. muons with opposite sign to the 
muons stored in the neutrino factory. Wrong-sign muons (fJ-~) result from v e — » oscillations (if /i + are stored), 
and can be used to measure the mixing angle #13, determine the neutrino mass hierarchy, and search for CP violation 
in the neutrino sector. In addition to the wrong-sign muon signal, there will also be right-sign muon events. These 
events come from the disappearance muon neutrino channel, (y^ — > u^), if positive (negative) muons are 

stored. The discrimination of the wrong and right sign muons requires the identification of charged current (CC) 
muon neutrino interactions, and the measurement of the sign of the produced muon. If the interacting neutrinos 
have energies of more than a few GeV, standard neutrino detector technology, based on large magnetised sampling 
calorimeters, can be used to measure wrong-sign muons with high efficiency and very low backgrounds. This has been 
shown to work for neutrino factories with energies of about 20 GeV or greater 1(| ID, [28| . 



Lower energy neutrino factories [2_4j,|27|, which store muons with energies < 10 GeV, exploit a fully active calorimeter 
within a magnet, a detector technology which ensures the detection of lower energy muons. The possibility of a low 
energy neutrino factory with non-magnetic detectors has also been explored in [29( , although we have not considered 
it here. A neutrino factory with muon e nerg ies of about 4 GeV has been shown to enable very precise measurements 
of the neutrino mixing parameters [13, Hjf ■ Electron charge identification also becomespossible in a low energy 
neutrino factory equipped with a magnetised totally active scintillating detector (TASD) [28j . Therefore, in addition 
to the wrong and right-sign muons, there will also be wrong and right-sign electrons from the appearance channel (the 
platinum channel), — > v e , and the disappearance channel, v e — > v e , for positive muons stored in the decay ring . 
These platinum channels, which are the T-conjugates of the golden channels, could provide a possible way of resolving 
the problem of degenerate solutions [30M33I ] . It is well known that even a very precise measurement of the appearance 
probability for neutrinos and antineutrinos at a fixed L/ E allows for different solutions of (#13, sign(Am| 1 ), 8), severely 
weakening the sensitivity to the se p arameters. Many strategies have been advocated to resolve this issue which in 
general involve another detector [34M39j , the combination with another experiment (22], [23], l40l - |49l | and / or the addition 
of new channels [ill H3| ■ 

We will consider the impact of the addition of the platinum channels in both a low statistics and high statistics 
scenario. We will show that if the platinum channel is accessible with negligible backgrounds (by which we mean a 
charge misidentification rate of < 10~ 2 , and that fewer than 10 -2 of all neutral current events are wrongly counted 
as signal events - see Section Hi Bl for more details), in the case of low statistics, its addition would greatly increase 
the potential of the low energy neutrino factory. However for the high statistics scenario, the platinum channel has 
much less of an effect - higher statistics combined with complementary information from the different energies of a 
broad beam alone are sufficient to resolve degeneracies and maximise the performance of the setup. Unfortunately, 
once realistic backgrounds from misidentified pions and the challenging electron charge identification are taken into 
account, the addition of the platinum channel does not provide any significant improvement. 

The structure of the paper is as follows: in Section |H] we discuss in detail the physics reach of the proposed setup, 
which exploits the wrong and right sign muon and electron signals. Based on the reasons described in Section [Hi the 
decaying muon statistics assumed here is higher than in previous studies. The assumptions for the detection efficiencies 
and energy resolution of the detector have also been modestly improved, based on the simulated performance of 
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the NO^A TASD detector |7|. These assumptions provide a more competitive setup with respect to our previous 
studies [24], [27} ■ We perform detailed numerical simulations and discuss the sensitivity of the low energy neutrino 
factory to the mixing angle #13, to the CP violating phase 6, to the neutrino mass hierarchy, to the octant of the 
atmospheric mixing angle #23 and to deviations from maximal atmospheric mixing as a function of the energy resolution 
of the detector and the number of muon decays per year (with and without the addition of the platinum channels). 
In Section Hm we introduce our preliminary studies of a magnetised 100 kton liquid argon (LAr) detector, comparing 
its performance to that of the TASD and other near term and future long-baseline neutrino facilities. Finally, in 
Section HVl we draw our conclusions. 

II. PHYSICS REACH: OPTIMISATION OF THE EXPERIMENTAL SETUP 

In this section we present the results from numerical simulations of the low energy neutrino factory. We have used 
the GLoBES software package [5l| to simulate several experimental configurations which will be described in the 
following subsections. These have led us to an optimised setup, which we use unless otherwise specified, defined by 
the following: the baseline is 1300 km, corresponding to the Fermilab to DUSEL distance. For the beam we consider 
a muon energy of 4.5 GeV with 1.4 x 10 21 useful muon decays per year per polarity, running for ten years. This flux 
is larger than the 2.5 x 10 20 useful muon decays per year per polarity usually considered for the International Design 
Study (IDS) neutrino factory [28] for three reasons. Firstly, the IDS neutrino factory distributes the beam to two 
baselines whilst only one baseline is required for the low energy neutrino factory - therefore a factor of two is gained. 
Moreover, 1 x 10 7 operational seconds per year were assumed in 28] whilst we consider that 2 x 10 7 operational 
seconds per year should be achievable. However, in order to perform an equal comparison with other facilities, 1 x 10 
seconds per year were assumed in the comparison plots of Section UTTl Finally, the extra factor of 1.4 arises from a 
re-optimisation of the accelerator complex for the lower muon energy required with respect to the IDS design [52j . 

For the detector we assume a totally active scintillating detector (TASD) with a fiducial mass of 20 kton, energy 
threshold of 0.5 GeV, energy resolution of 10% with 19 variable-width bins, efficiency for /i ± detection of 73% below 
1 GeV and 94% above, and a background level of 10~ 3 on the v e — > (p e — > v^) and — > (y^ — > D^) channels. 
We assume that the background to each channel arises predominantly from charge misidentification and neutral 
current events, modeling the background to each channel as a constant fraction of the rates of the wrong-sign and 
neutral current channels. For systematics, we use 2% on both the signal and background, assuming the errors to be 
uncorrelated. 

For the platinum (y^ — > v e and — ¥ D e ) channels, identifying the charges of the electrons at the energies involved 
will be very challenging. In addition, the pion background is very difficult to separate from the electron samples 
as compared to the muon samples. Preliminary estimations [53[ suggest that the efficiency for e ± detection could 
be 37% below 1 GeV and 47% above, and that the backgrounds (we make similar assumptions for the background 
sources as for the muon channels) could be reduced to a level of 10~ 2 . Although it is uncertain as to whether this 
detector performance can actually be achieved - further simulations are required - we will regard these numbers as an 
optimistic estimate for the electron detection capabilities of a TASD. 

We assume the same oscillation parameters as in sin 2 6*12 = 0.3, #23 = 7r/4, Am 2 ! = 8.0 x 10~ 5 eV 2 , and 

IAto^I = 2.5 x 10 -3 eV 2 with a 10% uncertainty on the atmospheric parameters, 4% uncertainty on the solar 
parameters, and 2% uncertainty on the matter density. In all our simulations we have used the exact oscillation 
probabilities, taking into account matter effects, and have marginalised over all parameters. 

The most significant alteration relative to the previous setup 0, H3] consists of the addition of the — > v e and 

— > D e channels. These platinum channels are the T-conjugates of the golden channels (y e — > Uu, D e — > P M ), and we 
will investigate the power of this combination to reduce the degeneracies in the 6>i 3 , 5 and sign(Am 2 1 ) parameter space. 
It has been shown that typically, the elimination of these degenerate solutions may require additional information 
from a second baseline and detector [ijj |3~4T - [39T | , or from a complementary experiment [13, HH, [40l - |49l | . 

We try instead to exploit the ability of the TASD to detect and identify the charge of e~ and e + , which gives access to 
the platinum channel. The probability for this channel, to leading order in the small quantities #13, a — Am^/Am 2 ^ 
and EA/Am^ (where A = \J2G-p'n e is the matter potential and n e is the electron number density), is identical to 
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that for the golden channel [16[ with the interchange of S — >• —5 and is shown below 



P»e = ^ 13 S2 3 sin 2 (^f -) (2-la) 
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1/'.' 2 
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We use a notation where = sin 6^, s 2i j = sin(20jj), Cjj = cos(9jj, c 2 ij — cos{29 i j) 1 E is the neutrino energy and 
i is the baseline. The first line of each probability, subequations (|2.1a[) and (I2.2a[) , is the atmospheric term which is 
quadratic in sin(26*i3) and will be dominant in the scenario that 9\ 3 is large (sin 2 (26'i3) > 10~ 2 ), and at high energies. 
The atmospheric term provides sensitivity to #13, the mass hierarchy, and is sensitive to the octant of #23- The second 
line, subequations (|2.1b[) and (|2.2b[) . is the CP term which is linear in sin(26*i3) and dominates for intermediate values 
of 0i3 if 5 7^ or 7r. The dependence on S enters via the oscillatory cosine term which can take either a positive or 
negative sign depending on the value of the phase. This can lead to constructive or destructive interference between 
the atmospheric and CP terms, meaning that sensitivities to #13 and the mass hierarchy are strongly dependent 
on the value of 5. Due to the inverse dependence on energy, the CP term becomes most visible at lower energies; 
therefore it is important to have access to the second oscillation maximum to establish if CP is violated. Thus a 
lower energy is desirable to enable a clean measurement of 5, whereas a higher energy and, especially, a long-baseline, 
guarantees sensitivity to the mass hierarchy. The low energy neutrino factory is unique in having a surprising degree 
of sensitivity to the mass hierarchy in spite of its low energy (as we show in Section III CP due to its broad spectrum 
that includes energies beyond the first oscillation peak, thus enabling complementary information to be obtained to 
resolve degeneracies. The third line, subequations (I2.1cl) and (|2.2c|) . is the solar term which is independent of #13, 
S and the mass hierarchy, and is dominant in the case that #13 is very small (sin 2 (26*13) ^5 10~ 4 ). In this re gim e, 



measurements will be extremely challenging and a high energy neutrino factory may be the only option [16j, [18, 28]. 

If we consider the fact that the probability for the CP-conjugated golden channel, D e — > D^, takes a similar form 
to that of the golden channel but with the substitutions 5 — > — 5 and A — > — A and that the CPT-conjugated golden 
channel is identical to the golden channel, with the exchange of A — >• —A, we can understand the complementarity of 
these four channels: each of the channels has a different dependence on the parameters 6*13, 5 and sign(Am§ 1 ) and 
so degenerate solutions are present at different points in the parameter space for each of the channels [50j . Thus the 
degenerate solutions from one channel can be eliminated by the information from another channel. 

We will mention briefly that the ability of the TASD to detect electrons also enables measurement of the v e (D e ) 
disappearance channel whose probability is 

9 o / Amli AL\ 
P„._>„. « l-s^sm 2 -J . (2.3) 

However, as this channel is CP-invariant and has only a weak dependence on the mass hierarchy, it is expected and 
has been verified in our study that its addition does not provide any significant improvement. 

In the rest of this section we now show the impact of our improved statistics and energy resolution, and of the 
addition of the platinum channels. 

A. Energy resolution 



Firstly we illustrate how our more optimistic estimate of 10% for the energy resolution, dE/E, (with 19 variable 
width bins) improves the performance of our setup. In Fig. Q] we show how the new resolution improves upon the 
old value of 30% (with nine variable width bins) for 6*13 = 1° and 5°. We observe that in addition to the significant 
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dE/E = 30% 
dE/E = 10% 



10" 4 10" 3 10" 2 10" 2 10" 1 10° 

sin 2 (29 13 ) sin 2 (29 13 ) 

(a) Sis = 1° (b) 013 = 5° 

FIG. 1: Comparing an energy resolution of dE/E = 30% (dotted blue lines) and 10% (solid red lines): 68%, 90% and 95% 
confidence level contours in the sin 2 (2#i3) — 8 plane for true values of 8 = —180°, —90°, 0° and 90° and a) #13 = 1°, b) #13 = 5°. 

increase in sensitivity to #13 and 5, the hierarchy degeneracy is now almost completely eliminated even for small values 
of 813. This is possible because at these low energies, the oscillation probability displays a rich dependence on the 
oscillation parameters as a function of the neutrino energy. In particular, the hierarchy can be determined from the 
position of the first oscillation maximum, and the value of S from both the first and second oscillation maxima. Hence 
if the neutrino energy can be measured with sufficient precision to enable the oscillation spectrum to be accurately 
reconstructed, a significant improvement in the sensitivities to the oscillation parameters can be achieved. 

B. Inclusion of the platinum channel 

We define Scenario 1 to be the one in which only detection is possible, giving us access to only the and 
appearance and disappearance channels. In Scenario 2 it is also possible to detect e ± and hence exploit the additional 
information from the v e and v e appearance channels. To illustrate the impact of the addition of the platinum channel 
to our setup, in Fig. [5] we compare the sensitivities of the two scenarios when using a muon decay rate of 5.0 x 10 20 (left 
column) and 1.4 x 10 21 (right column) per year, varying the background level of the v e (i> e ) appearance channel from a 
hypothetical zero (top row) to 10~ 2 (bottom row). By this we mean that we have simulated a charge misidentification 
rate such that 10 -2 of all e~ (e + ) are wrongly identified as e + (e~), and that 1CP 2 of all v e (p e ) neutral current events 
are wrongly counted as signal events. 

In the case of the lower statistics, we observe that the addition of the platinum channel with zero background 
produces a drastic improvement in sensitivity to all parameters. For a background of 10 -2 the improvement is much 
smaller but can still help to alleviate the hierarchy degeneracy (see [H(J). At higher backgrounds we find that this 
gain is lost. In the case of the high statistics, we already observe a smaller improvement for zero background, which 
becomes insignificant at a background level of 10~ 2 . Thus we conclude that since the estimated background on the 
v e (z/ e ) appearance channels will be at best ~ 10~ 2 , the platinum channel could help in the measurement of the mass 
hierarchy if statistics are limited to 5.0 x 10 20 useful muon decays per year, whereas it will be almost irrelevant for the 
higher statistics scenario. An increase in statistics in the golden channel thus provides a much larger improvement 
in the facility performance than the addition of the platinum channel if background levels below 10 ~ 2 cannot be 
achieved. 



C. Sensitivity to #13, CP violation and mass ordering 

Here we present the results of our optimisation studies, in terms of 3cr $13 discovery potential, CP discovery 
potential, and sensitivity to the mass hierarchy in the sin (2^13) — 8 plane (Fig. [3]). In addition we also consider the 
3cr sensitivity to #23 in the sin 2 (2^13) — sin #23 plane, both in terms of the ability to exclude a maximal value of #23 
(Fig. 0k) an d to identify the octant of 6*23 (Fig. \Qp). The results from our optimised setup described in Section [TTI are 
shown by the solid green lines. We have also considered a setup where only the statistics are altered, to 2.8 x 10 21 



6 




sin 2 (29 13 ) 



(a) 5.0 X 10 20 decays per year, zero background 




-150 - 

10" 4 10" 3 10" 2 

sin 2 (29 13 ) 

(b) 1.4 X 10 21 decays per year, zero background 




sin 2 (29 13 ) sin 2 (29 13 ) 



(c) 5.0 X 10 20 decays per year, 10 2 background (d) 1.4 X 10 21 decays per year, 10 2 background 

FIG. 2: Comparison of Scenario 1 (u^ (v^) appearance and disappearance only - dotted blue lines), and Scenario 2 (y e (p e ) 
appearance included - solid red lines) when using 5.0 x 10 20 (J, decays per year (left) or 1.4 X 10 21 decays per year (right), and 
a background of zero (top row) or 10 -2 (bottom row) on the u e (u e ) channels: 68%, 90% and 95% confidence level contours in 
the sin 2 (26»i 3 ) - S plane, for S = -180°, -90°, 0° and 90° and 9 13 = 1°. 

decays per year (solid red lines), and a setup where only the muon energy is increased to 6.0 GeV (dashed blue lines). 
From this we demonstrate that for all the observables considered, doubling the flux is always preferable to an increase 
in energy. 

For $i3 discovery potential, CP discovery potential and 623 sensitivity we only show the results for a normal 
hierarchy, having verified that similar results are obtained for an inverted hierarchy. We have assumed in Fig. |4] 
(6*23 sensitivity) a value of S = 90° although we have also studied other values of 5 and find no strong dependence 
on the CP phase, since sensitivity to #23 is mainly obtained from terms with no dependence on S in the oscillation 
probabilities discussed in Section [TTJ For the exclusion of #23 = 45°, an upward curve is seen for large 813. This can 
be understood because the addition of a large 613 to the disappearance probability introduces an asymmetry in 
$23 that shifts the contours to larger values (see eq. (1) and Fig. 8 of Ref. |55j|). 

We note that this setup has remarkable sensitivity to 6*13 and S for values of sin 2 (2#i3) > 10~ 4 , and that its sensitivity 
to the mass hierarchy is an order of magnitude better that that of other proposed experiments exploiting the same 
baseline e.g. the wide-band beam experiment in {56[. We can attribute these qualities to the unique combination of 
high statistics and good background rejection coupled with an intermediate baseline, allowing for a clean measurement 
of the CP phase whilst also allowing for the mass hierarchy to be determined for sin 2 (20i3) > 1CP 3 . 

We have also explored how the precision with which #13, S and the deviation from maximal #23 could eventually be 
measured at the low energy neutrino factory, varies as a function of exposure (detector mass x decays) per polarity. 
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FIG. 3: 3a confidence level contours in the sin 2 (2#i3) — S plane for a) #13 discovery potential, b) CP discovery potential, c) 
hierarchy sensitivity (for true normal hierarchy), d) hierarchy sensitivity (for true inverted hierarchy). 



Our standard setup corresponds to 20 kton x 1.4 x 10 21 decays/ year x 10 years = 2.8 x 10 23 kton x decays per 
polarity. We find that if the mixing angle #13 turns out to be large, the unknown leptonic mixing parameters could be 
measured with unprecedented precision at a future low energy neutrino factory for sufficiently high exposures. The 
gain in precision is much less pronounced for values larger than 6 x 10 23 kton x decays per polarity, hence it may not 
be worth trying to increase the exposure beyond this value. 

Fig. [5ji shows the la error expected in the measurement of the mixing angle #13 at a future low energy neutrino 
factory as a function of the exposure (in kton x decays) per polarity, assuming that nature has chosen #13 = 5°. The 
dependence of these results on the value of the CP violating phase is very mild. The la error in the extraction of #13 
when no backgrounds and no systematic errors are included in the analysis is illustrated in Fig. [SJd. Comparing the 
two panels we observe that non-zero systematics and backgrounds effectively halve the exposure. 

Fig. [5^, shows the la error expected in the measurement of the CP phase S as a function of the exposure for a 
simulated value of #13 = 5°, for different values of the CP violating phase 5. The results are highly dependent on the 
value of the CP violating phase, as expected. For S = 90°, there are strong correlations with 6*13, as can be seen from 
Fig. [21 and therefore the error in the measurement of the CP violating phase is larger. The la error in the extraction 
of 5 when no backgrounds and no systematic errors are included in the analysis is illustrated in Fig. |5}d. Switching 
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FIG. 4: 3a allowed regions in the sin 2 (2^13) — sin #23 plane for a) potential to exclude #23 = 45°, b) sensitivity to the #23 octant. 
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FIG. 5: lcr error in the measurement of the 613 mixing angle for a simulated value of 813 — 5° and different values of the CP 
violating phase S when a) including systematic errors and backgrounds, b) no systematic errors and backgrounds are included. 



off systematic errors and backgrounds has a larger impact for the 6 = 0° case, again effectively halving the exposure, 
since correlations among 5 and #13 are negligible when 6 = 0° and the precision is more limited by the background 
and systematic errors instead. 

We also explore the sensitivity to maximal mixing, i.e. the ability to exclude #23 = 45°, versus the exposure. We 
present the 3a results in Fig. [7] We have used a simulated value of 6*13 = 0° here (so that 5 is irrelevant) as the 
sensitivity to 623 maximality comes primarily from the (p^) disappearance channels which are not dependent on 
6*13. Since the disappearance channels are also not strongly dependent upon systematic errors or backgrounds, there 
is little change obtained by switching these off. 



9 




- 6=0° 

- 5=90° 
90 



2.X10 23 4.xl0 23 6.X 10 23 8.X 10 23 1.x 10 24 
kton-decays 

(a) Sensitivity to 8, including systematic errors and 
backgrounds 



2.X 10 23 4.X 10 23 6.xl0 23 8.X 10 23 1.x 10 24 

kton-decays 

(b) Sensitivity to <5, no systematic errors and 
backgrounds 



FIG. 6: la error in the measurement of the CP violating phase <5 for a simulated value of #13 = 5° and different values of 
the CP violating phase S when a) including systematic errors and backgrounds, b) no systematic errors and backgrounds are 
included. 
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FIG. 7: 3a regions for which maximal 823 can be excluded, using a simulated value of 813 = 0° when a) including systematic 
errors and backgrounds, b) no systematic errors and backgrounds are included. 



III. LIQUID ARGON DETECTOR AND COMPARISON WITH OTHER EXPERIMENTS 

Recently there has been much interest in the possibility of constructing a kton-scale liquid argon (LAr) detector [{SlJ ■ 
If such a detector can be magnetised, it could be utilised in combination with a low energy neutrino factory and we 
have performed some preliminary studies to assess the potential of a 100 kton LAr detector for this experiment. As the 
design of large LAr detectors is still in the early stages, there are large uncertainties in the estimates for the detector 
performance. We assume an efficiency of 80% on all channels and 5% energy resolution for quasi-elastic events, then 
consider a range of values for other parameters. In the most conservative scenario, we assume 5% systematics, 20% 
energy resolution for non quasi-elastic events, and backgrounds of 5 x 10 -3 on the (vp) (dis) appearance channels 
and 0.8 on the v e (P e ) appearance channels [58|. For the optimistic scenario we use values identical to the TASD: 
2% systematics, 10% energy resolution for non quasi-elastic events, and backgrounds of 1 x 10~ 3 on the (v^) 
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(dis) appearance channels and 1 x 10 -2 on the v e (v e ) appearance channels. We find that varying the systematics, 
energy resolution and v e (v e ) background do not play a large role in altering the results; the dominant effect comes 
from the variation of the (v^) background. 

In Fig. [5] the results of the low energy neutrino factory with both the TASD and the two assumptions for the LAr 
detector are compared with other long-baseline experiments. We show the 3tr results for #13 discovery, CP discovery 
potential, and hierarchy sensitivity (for normal hierarchy only) as a function of sin 2 (2#i3) in terms of the CP fraction. 
In order to make a fair comparison, we have used half the flux described in Section [TT] for the low energy neutrino 
factory, to make it consistent with the other experiments which assume only 10 7 seconds per year of observation. 
However, we believe that the fluxes in Section [TT] are feasible. The results from the TASD are shown by the red line 
and those from the LAr detector are shown by the blue band. The right-hand edge of the band corresponds to the 
conservative estimate of the detector performance, and the left-hand edge to the most optimistic estimate. As the 
optimistic scenario assumes an almost identical performance to the TASD, the left-hand edge of the blue band also 
corresponds to the results obtainable from a 100 kton TASD. Results from the high energy neutrino factory [28j j . 
wide- band beam [56j], T2HK (59j | . IOO7 /3-beam |60j . 3507 /3-beam [6l| and 4- ion /3-beam [62| are also shown. 

In terms of sensitivity to #13, a conservative low energy neutrino factory is an order of magnitude less sensitive 
than the high energy neutrino factory, but is still competitive with the /3-beam experiments, giving an approximately 
equal performance to the 4-ion /3-beam (which requires two baselines to resolve the degeneracy problem, as for the 
high energy neutrino factory). However, the performance of an aggressive low energy neutrino factory setup surpasses 
that of all other experiments except for the high energy neutrino factory. For CP violation, the low energy neutrino 
factory gives remarkable results: the most optimistic setup outperforms the high energy neutrino factory for all values 
of 0i3, and even the most conservative setup gives a superior performance for sin 2 (20i3) > 2 x 10~ 3 , again equaling 
that of the 4-ion /3-beam. For sensitivity to the mass hierarchy, the low energy neutrino factory gives an improvement 
over all other experiments apart from the higher energy setup and the 4-ion /3-beam with their challenging 7000 km 
baseline. 



IV. CONCLUSIONS 

We have optimised a low energy neutrino factory setup with a baseline of 1300 km, defining a reference setup to be 
one with a muon energy of 4.5 GeV and 1.4 x 10 21 useful muon decays per year, per polarity, running for ten years. For 
the detector we assume a totally active scintillating detector (TASD) with a fiducial mass of 20 kton, energy threshold 
of 0.5 GeV, energy resolution of 10%, efficiency for ^ detection of 73% below 1 GeV and 94% above, efficiency for 
detection of 37% below 1 GeV and 47% above, and a background level of 10~ 3 on the v e — f (p e — > P M ) and 
v yt Vfi SyO channels and 10 -2 on the — > v e (v^ — > u e ) channels. We have also considered a 100 kton liquid 

argon detector and found that its performance would equal or surpass that of the 20 kton TASD. 

We have demonstrated how improving the energy resolution and statistics improves the performance of the setup, 
showing that in particular high statistics play a vital role. We have also shown how the combination of golden 
and platinum channels could be a powerful way of resolving degeneracies, especially in the case of limited statistics. 
However, once realistic background levels of at least 10 -2 are considered, the improvement achieved by adding the 
platinum channel is negligible. Therefore, more effort should be invested into achieving larger statistics for the golden 
channel than in improving the platinum channel, at least for standard physics searches. 

Using our optimised setup, the low energy neutrino factory can have sensitivity to #13 and 5 for sin 2 (2813) > 10~ 4 , 
competitive with the high neutrino factory. Sensitivity to the mass hierarchy is accessible for sin 2 (2#i3) > 10 -3 , 
better than other experiments using the same baseline due to the complementarity of measurements with different 
channels and different energies. Even if the flux is halved to equal that of other long-baseline experiments, the low 
energy neutrino factory is still competitive, performing especially well for CP discovery at large values of #13. We 
have also studied the sensitivity to 823, finding that it is possible to exclude maximal 823 at 3u for 823 ^ 43° and 
823 ^ 47°, roughly independent of #13, and to identify the octant for 823 ^5 37° and 823 ^ 53°. 

Studies of the sensitivities as a function of exposure (detector mass x number of decays) show that the effect of 
non-zero systematic errors and backgrounds is to effectively halve the exposure, affecting the sensitivity to #13, S 
(especially for S = 0°) and 823- For exposures > 6 x 10 23 kton x decays per polarity and large 813, the low energy 
neutrino factory could measure the oscillation parameters with unprecedented precision. 

We conclude that the low energy neutrino factory has excellent sensitivity to the standard oscillation parameters 
and is therefore a potential candidate for a next-generation long-baseline experiment. 
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(c) Hierarchy sensitivity 

FIG. 8: Comparison of 3u allowed contours for the low energy neutrino factory with a 20 kton TASD (red line) and 100 kton 
LAr detector (blue band), the high energy neutrino factory (black line), the wide-band beam (purple line), T2HK (yellow 
line) and three /3-beams (green, orange, light blue lines) for a) 6*13 discovery potential, b) CP discovery potential, c) hierarchy 
sensitivity. 
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